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Abstract 

We consider the four dimensional Euclidean Maxwell theory with a Chern-Simons term on 
the boundary. The corresponding gauge invariant boundary conditions become dependent on 
tangential derivatives. Taking the four-sphere as a particular example, we calculate explicitly 
a number of the first heat kernel coefficients and obtain the general formulas that yields any 
desired coefficient. A remarkable observation is that the coefficient 02, which defines the one- 
loop counterterm and the conformal anomaly, does not depend on the Chern-Simons coupling 
constant, while the heat kernel itself becomes singular at a certain (critical) value of the coupling. 
This could be a reflection of a general property of Chern-Simons theories. 
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1 Introduction 



The heat kernel expansion plays an important role in quantum field theory. In fact, the ultraviolet 
divergences and renormalization structure of the theory at the one-loop approximation are 
encoded in a few of the first heat kernel coefficients a n . For massless theories in d dimensions the 
only divergent term is proportional to 0,4/2- The heat kernel expansion is closely related to the 
quantum anomalies, most notably, to the conformal and chiral anomalies. It has also numerous 
applications in mathematics. No wonder that much effort has been devoted to calculate the heat 
kernel expansion. For operators of the Laplace type, the coefficients a n are known fairly well, 
both for manifolds without a boundary and for Dirichlet and Neumann boundary conditions in 
manifolds with a boundary. More details on the physical applications of the heat kernel expansion 
can be found in the recent monographs Q. The mathematical background is described in ||. 

During the last couple of years, considerable interest has been attracted by the heat kernel 
expansion for boundary conditions depending on tangential derivatives @-f|, following the 
pioneering work M. In the mathematical literature such boundary conditions are called oblique. 
For physicists, the interest towards oblique boundary conditions has been primarily motivated 
by problems of quantum gravity |10|| . It has been realized that the diffeomorphism invariant 



boundary conditions for the graviton inevitably contain tangential derivatives [11, H, 12, 13 



Oblique boundary conditions may also appear in electrodynamics, when the properties of a 
physical boundary depend on the photon frequency. A more formal motivation for studying 
oblique boundary conditions can be found in fTJJ, where tangential derivatives arise when one 
calculates vacuum expectation values of second order differential operators. 

The usual strategy @]-|§ employed in the calculation of the heat kernel coefficients for 
oblique boundary conditions is as follows. By means of the corresponding invariant theory, each 
coefficient can be expressed through several universal functions, which are calculated by using 
ordinary conformal variation techniques and some explicit examples. However, up to now, only 
a very few examples admitting calculation of the spectrum and an explicit evaluation of the 
heat kernel were known || ||, In the present paper we study the Euclidean Maxwell theory 
with a Chern-Simons boundary term, which generates oblique boundary conditions. We find 
the spectrum of the Laplace operator on a ball and calculate the heat kernel expansion. 

Apart from its application to the theory of the heat kernel expansion, our model may be 
interesting by itself. Tangential derivatives appear on the boundary conditions in a very natural 
way. An intriguing property of the model, established in the present paper, is that the conformal 
anomaly and the one-loop counterterm do not depend on the charge standing in front of the 
Chern-Simons action. Other heat kernel coefficients become singular for a certain value of this 
charge, what suggests the existence of a critical value of the Chern-Simons charge. Let us recall 
that the Chern-Simons gauge theory introduced in \TE\ exhibits some unusual and intriguing 
properties, as, for example, generation of states with fractional statistics first observed by Wilczek 
p| . Later on this phenomenon was used in the theory of the Fractional Quantum Hall Effect 



(FQHE) (for a recent review see JIT)). In spite of the fact that the geometry of our model is 
different from the ones usually considered in the FQHE, the results of our calculation may be 
important in order to gain insight into the general properties of Chern-Simons theories, as e.g. 



non-renormalization theorems 118 



The paper is organized as follows. In section 2 we formulate the model, discuss its general 
properties and define the eigenf unctions on the Euclidean ball. In section 3 we calculate the heat 
kernel asymptotics explicitly. A discussion on the possible significance of the results obtained (a 



2 



list of the first ten heat kernel coefficients is given in Table 1) is presented in section 4. 

2 Chern— Simons boundary conditions 

Consider the action for the Euclidean Maxwell theory on a manifold M 

S= \ I d A x g^F^F^, F^ = d^A v - d v A„. (1) 

We can add to the action ([I]) a boundary term. If we require gauge and coordinate invariance to 
hold and, if we do not want to introduce any dimensional parameter, the only choice available 
is the Chern-Simons action 

Scs = ^l d 3 x e^AdjAk, (2) 

2 JdM 

where e l ^ k is the Levi-Civita tensor and the x\ j = 1, 2, 3, are coordinates on the boundary dM, 
a being a real parameter. 

To calculate the path integral, it is convenient to write S + Scs i n the form / ALA, with L 
a second order differential operator. To this end, we integrate by parts, obtaining 

S + Scs = ~/ ^jsV-rA + VV^ 

2 JM 

+1 f ^(h^dtfAi - diA N )Ai + ae^AdjAk), (3) 
2 JdM 

where A is the Laplace operator, N the outward pointing normal vector, and h is the determinant 
of the induced metric on dM. There are, at least, two sets of gauge invariant boundary conditions 
which ensure vanishing of the surface term in (^|). The first set is ordinary relative (or magnetic) 
boundary conditions: A^qm = 0, i = 1, 2, 3, (d]y + k)A N \Q M = 0. Here k is the trace of the second 
fundamental form of the boundary. These boundary conditions have been extensively studied 
in the literature, but we shall not consider them here. Another possible set is the following: 

A N \ 9M = 0, {d N A i + ah-hj k d j A k )\ 9M = 1 = 1,2,3. (4) 

For a = 0, Eqs. (|]) become ordinary absolute (or electric) boundary conditions. 

The boundary conditions (£|) possess two properties which make them useful for quantum 
electrodynamics. First, they are gauge invariant. This means that if A^ satisfies (Q), then A + <90 
also satisfies (|), provided that 

d N 4>\dM = 0. (5) 

After quantization, Eq. (|) becomes the boundary condition for the ghost field. Second, the 
Laplace operator is symmetric, i.e., 

/ d A x gtfAWfAA® - A^^AA^) = 0, (6) 

JM ^ * 

if both A^ and A^ satisfy (|). 

It is natural to impose the Lorentz gauge condition 

VM„ = 0. (7) 
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In this case the path integral is given by 

Z = det~'(-A) x detf(-A), (8) 

where the first determinant is calculated on the space of transversal vectors, with the boundary 
condition (f|), and the second one on the space of scalar fields satisfying (H). 
As an example, consider a ball with the metric 

ds 2 = (dx ) 2 + (x°) W, < x° < r, (9) 

where dfl 2 is the metric on the unit sphere S 3 . 

We can use the basis of Ref. |19j in the space of transversal vector fields: {A T } = {A L , A(ip)}, 



where 

Aq = 0, (3) V^ = 0, 

Aoty) = (3) Ar^, A J (^) = - (3) V l (9 + -)r^, 2 = 1,2,3, (10) 

r 

ip being a scalar field, and and ®A the covariant derivative and the Laplacian on S 3 , 

respectively. 

The boundary conditions for the field i[) do not depend on a. Hence, the contribution of ip to 
the heat kernel and functional determinant is the same as for the absolute boundary condition 
[ [B]| , EDI. Therefore, let us concentrate on the A 1 - contribution. The operator h~^£^ k dj can 
be diagonalized on the unit sphere S 3 (see, e.g., [plf ). It has for eigenvalues ±(Z + 1), with 
degeneracies given by 

Df = 1(1 + 2), 1 = 1,2,.... (11) 

The corresponding eigenvalues of the Hodge-de Rham Laplacian on S 3 are — (/ + l) 2 . The 
eigenf unctions of the vector Laplace operator on the unit ball are found to be 

J m (Ax°)r / ± (x i ), (12) 

where Ji + i are ordinary Bessel functions and Y l are vector spherical harmonics. The eigenvalues 
of the Laplacian are —A 2 , where the A's are defined, implicitly, through the equation for the 
boundary condition: 

XfJ f l+1 (Xt) ± J I+ i(Af ) = 0. (13) 

r 

The prime denotes here differentiation with respect to the argument. The contribution of A 1 - to 
the path integral reads 

If one proceeds, as usually, by splitting the calculation of the determinant into two parts, one 
has to take into account a possibly non- vanishing determinant anomaly (see, e.g., 



3 The zeta function: calculation of the singularities 

In the following, we will use the connection between the zeta function and the heat kernel: 

oo 

CHs) = f f 7) Jdtt^j:e-^ )2 \ (15) 
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in order to obtain the heat kernel coefficients Cf or at'- 



(4tt) 2 
1 



(4V 



2 a n/2 



( 4vr )3/2 



E* 

c, 



n=0 
± 



U n/2 



(n-l)/2 



Res 



(16) 



The shifted heat kernel coefficients Cj°' = 



C~ are commonly used in Casimir energy cal- 

..... j_ " " 

culations. Calling the functions that determine implicitly the spectral values corresponding 



to the given boundary conditions 



the zeta function is obtained as 



oo oo „ 

c ± (^) = EA ± (Af)- 2s = E^ + 2) / 

1=0 1=0 T 



/ 7 2ttz dk 1 y 1 



(17) 



with 



(18) 



where we have transformed the sum over the spectrum into a contour integral which goes around 
all the real zeros of the function $ defining the boundary conditions. 

By doing this we have reduced the problem to the evaluation of certain contour integrals, 
which can be solved in an explicit — though technically quite involved — way. The interested 
reader is addresed to the Refs. p3[ p4p , where this powerful method is described in great detail 
— in the particular case a = 0. Having understood this case, however, the situation here is not 
difficult to grasp, the main differences being in Eqs. (^3|) and (|24"D to follow, which will lead, in 
the end, to more involved hypergeometric functions. In any case, the final expressions will be 
still explicit, allowing for the straightforward use of very quick algebraic computation machinery, 
to obtain the final expressions of the heat kernel coefficients to any order. This will demonstrate, 
once more, the power (and adaptability) of the zeta-function procedure to perform these sort of 
calculations (which are extremely cumbersome, by any means). 

By deforming the integration contour around the imaginary axis, in the usual way [23, 2\ 
we obtain 



sin(7rs) 



7T 



x roo 

V/(/ + 2) / dkk~ 
i=o Jo 



9 i 

X dk lU 



fc- (m) (l\ +1 (kr)±^±^I l+1 (kr) 



(19) 



Introducing into the equation the asymptotic expansions corresponding to / and /' when k, I 
oo, with z — kr/(l + 1) fixed, we can write the zeta function under the form 



c±(s) = si^ £ )g ;i(i + 2)( '<+I' 



-2s 



7T 



Z=0 



dzz 2s — {{I + 1) In \ - — 

o dz P ; I zt V 1 + 1 



+ In 



1 + E^UaJl + E ^ 



t (i + iy 



t v + !) j 



(20) 
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where the functions Uj(t) and Vj(t) give the asymptotics of the Bessel functions in the way 



l v {yz) 



I'Jvz) ~ 



i + E 



u k (t) 



k=l 



V 



1 e^(l + z 2 )4 



7TZ/ 



1 | y Vk(t) 
fe=i ^ 



(21) 



with t = l/y/l + z 2 and 77 = V 1 + z 2 + In [2/(1 + yl + z 2 )}. The first few coefficients are given 
in 



22| , higher coefficients are easy to obtain using the recursions 

f 



Vk(t) 



l -t 2 (l - t 2 )u' k (t) + lJdr (1- 5r 2 K(r), 




u k {t) +t{t 2 - 1) 



-U k -i{t) + 



(22) 



(see p3| , p4}j). It is convenient to write the last logarithmic function in the expression of C ± ( s ) 
above as a series 



In 



u « + x ) 



n=l 



Hl±at) + J2 F n(t)(l + l)- n , (23) 
(1 ±at) n ' 



with pf n (t) a polynomial in t of degree 4n, which has the form 



3/t 



j=0 



(24) 



The use of this expansion can be based on the corresponding one for the case a = (see Refs. 

[M| ) and will be fully justified by the subsequent calculations. For any heat kernel coefficient 
only a limited number of terms in ( f2~3|) are needed. The integrations are readily performed using 
the basic result (Re s < 1): 



00 Q 

dzz~ 2s — 
az 



[1 ± at) n 



V(l-s) 



±a n 



T(s + 



T{n) 
¥) 



E(=Fa) 

ft=0 



r(n + /i)r(s + 



fc+/A 
2 J 



/i! r to) 



r(¥) 



r( s + 



i) 



r(l) 



n + 1 n fc + 1 

, — h 1, s H 

2 ' 2 2 



TT I \ n + I n k 
H [\—'2' S+ 2 





3 


{ 




1 


k 


2' 





3 fc + 1 



(25) 



where H is the generalized hypergeometric function. Observe that for the sum over the two 
signs of a the result simplifies to twice the second term on the r.h.s. of Eq. (|25|). 

However, one has to take into account that, as it stands, C ± ( s ) is a divergent quantity. 



We shall extract the singular parts of the zeta function according to the prescription in [23 
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24|, namely, by splitting the zeta function into a regular and a singular part and adding and 



subtracting a number, say M, of leading terms of the assymptotic expansion: 



M 



(±(s) = Z\s) + £ A± a (s) 



(26) 



The functions Af A (s) are found to be: 

r 2 T(s - 1/2) 



A n ,a( S ) 



V?r(s + 1) 



[C(2s-3)-C(2s-l)] 



[C(2s-2)-C(2s) 



,2s 



2T(s 



[C(2 S -2)-C(2,)]E(T 



a 



r(a + /i/2) 



m + h/2y 

i((2s + n-2)-((2s + n)} 
zt< ^, ^ h T(n + h)T(s + 



(27) 



j=0 1 W h=0 



j+n+h \ 
2 I 



where £(s) is the Riemann zeta function. Notice that all the series that appear here are in fact 
of the form corresponding to Eq. ( p5|) and thus give rise to generalized hypergeometric functions 
(the result simplifies, when adding the ± contributions). Let us study the poles of at half 
integer numbers, s = 3/2 — k, k G N. For even n = 2p the poles of the gamma functions at 
s = 3/2 — k on the r.h.s. of ( ^7|) do not contribute to A tot = A + + A~ due to the property 
z nj = (~^y z n,j- F° r °dd n = 2p — 1 the Riemann zeta functions on the r.h.s. have poles at 
integer values of s only. From those expressions, the following residua are obtained (fee N) : 



ResA to J(l/2) 

Res Ag(3/2) = - 
Resi4*(l/2) = - 

Res<%(3/2-n) 



2?r' 

^.3 



-3/2 



-1/2 



„3-2n 



3)i 



r(3/2-n) 

3n-l 



\ -v ^2n,2j \ -» 



oo „2h 



a 2h T(2n + 2h) T(j + h + 3/2) 
{2h)\T{n + j + h) 



v Z2n,2 j+ i y oF- x r(2n + 2/i-l) r(j + + 3/2) 



% r(2n) 



a-2n 



3)i 



(2/i- l)!r(/i + n + j) 



71 > 1, 



Res^, a (l/2-n) 



r(i/2-n) 

3n-l 



oo _2/i 



a 2/ T(2n + 2/i) T(j + h + l/2) 
(2h)\T(n + j + h) 



z *n,2 j+ i g o 2h ~ 1 r(2n + 2fe-l) r(j + h + 1/2) 



- r(2n) 



Res (3/2 -*) 



j=U - v- h=l 
('_l N )fc4 r 3-2fc 



(2h-l)W(h + n + j) 



n > 1, 



r(3/2 - fc) 



[C(2n - 2Jfc) - C(2n - 2A; + 2)] 
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k y^ 1 4n-i,2j ^yT* (-l)" + ^(2n + 2/i-2)!a 2fe 

f^ (2n-2)! ^ ( 2 / l )!(A;-n-j-/i-l)!r(n + i + /i-l/2) 

+ fc £ 2 4n-i, 2j+ i ""g- 1 (-l)"+' +fc (2n + 2/1-2)! a^- 1 



(2n-2)! ~ (2/i-l)!(A;-n-j-/i-l)!r(n + j + /i-l/2) 

n < — 1, A; < 4n, 

Res<_ M (3/2 - k) = ^ ' _ [C(2n - 2*) - C(2n - 2A; + 2)] (28 



x 



f 3 ^ 1 4n-i2j 3n ^ _1 (-l) n+J+fe (2n + 2/i-2)!a 



2/i 



\p (2n-2)\ ^ (2h)\(k-n-j-h-l)\T(n + j + h-l/2) 

3n-2 r + _ 3n-j-l w+i+ft/n , oh _ r>\| „2A-1 



v 4n-i, 2j+ i ^fT 1 (-l) n+J+fe (2n + 2^-2)!a 2 

(2n-2)! ^ (2/i-l)!(A:-n-i-/i-l)!r(n + i + /i-l/2) 

n < k — 1, k > An. 

Recall that the Riemann zeta function at the negative integers is given by Bernoulli numbers: 
C(—n) = —B n+ i/(n + 1). Using then the relation between these residua and the heat kernel 
coefficients 

fe-i 

ResC ± (3/2-A;) = Res J2 A L-i(3/ 2 - k ) + ResA^._ 2 (3/2 - k) + Res (3/2 - k) 

n=0 

(47r) 3 / 2 r(3/2 - k)' ^ 

we can calculate the C£ ot , for which an explicit, closed expression can be given, in terms of 
generalized hypergeometric functions of the type of Eq. ( p5| ) (what is quite clear from the 
consideration above). The formula can be immediately turned into a computer code to perform 
the final part of the calculation. It is important to observe that we have obtained, in fact, a 
closed, explicit formula for the heat-kernel coefficients, even if the actual calculation is much more 
involved than that for the case of the Laplacian on the sphere with Robin boundary conditions. 
With the help of a fast computer we can construct (in a couple of hours) a table of heat kernel 
coefficients Cl ot for integer values of k up to any (reasonable) order. The first few of them are 
given in Table 1. They are given as power series on a with rational coefficients. 

In a similar way we can proceed with the calculation of the coefficients for half-integer index 
Cfc+i/2' k £ ~N. The zeta functions C ± ( s ) have only one pole at integer s, with an a-dependent 
residue. In fact 

ResAf (1) = (30) 
cr — 1 

The values (±(1 - k) are needed and one has M = 2k. Here the Af(s) for i odd, % — 2j — 1, 
j G Nq, do not contribute. The relevant A\ ot (s) read now 



K\l-k) = r — [C{ _2k)-((2-2k)} 



2-2k 



K!aiX) = - T i 



r 2 / 7r 2 \ a 2 



2 3 J 1- a 



2 ' 



8 



A tot (1 k) r 2-2 k( lW , 1 v 3 f 1 f 4 fc+ i,2, f a 2fe (2fc + 2h)\T(j + h + 3/2) 

^ a 2fe - 1 (2A; + 2h - 1)! T(j + + 3/2) 1 
(2A;)! ^ (2/i- l)!r(fc + j + /i + l/2) J ' 

At* ri m r 2- 2fcf n * f . lV 3 y f 4-i, 2j f U 2h (2fc + 2ft-2)!rC7 + ft + l/2) 

- / (-1) (/c-lj. 2, i (2Jb_ 2 )!.^ (2/ l )!r(fc + j + / l -l/2) 



j=0 V. V'" h=0 



4-1^+1 g a 2fe - 1 (2fc + 2/i-3)!r(j + /i + l/2) | ^ 



(2ife-2)!^ (2/i-l)!r(fc + j + /i-l/2) 

And, for n G N, n < A; — 1, 

Ag(l-k) = -4r 2 - 2fe (fc-l)![C(2n-2A:)-C(2n-2A; + 2)] 



2h 



X 



min(fc-n-l,3n) ^± ^ min(A:-n-j-l,3 n-j) ^n+j+h^ + 2 h - 1)! a 

^ (2n-l)! ^ (2/i)!(A;-7i-j-/ l -l)!(7i + j + / l -l)! 

min( fe -n^2,3n-i) z + 2 _ +i mincfc-ng-LS n-j) (_i)n+ J +^(2n + 2fe - 1)! a 2 ^ 



(2n-l)! ^ (2/i-l)!(A;-n-i-/ i -l)!(n + i + / i -l)! 

Notice again that the summation ranges are different depending now on the fact that k < 4n + 1 
or k > An + 1. The corresponding heat-kernel coefficients are readily calculated from 

C^l - *) = £ ^(1 - *) + - k) + Af k+1 (l - k) = ( ~ 1)fc " 1( . fc ,~ 1)! C fc ± +1 . 

n=0 (47T)" 



2 



A number of the coefficients (7*^1 are listed in Table 1 . If desired, the coefficients a n can be 

then obtained by using Eq. fllEp. Let us finish this section by noting that here we have calculated 
the a-dependent part of the full heat kernel expansion for the Euclidean Maxwell theory with 
the boundary conditions (§) only Complete result can be immediately recovered by taking C tot 
from Table 1, subtracting the a part, and adding this to the known results [pi], [27], for 
absolute boundary conditions. All coefficients, including C±, are non-singular at a = 0. 



4 Conclusions 

The conditions @ we have considered in this paper are mixed and contain tangential deriva- 
tives. For boundary conditions of this type only the first non-trivial coefficient Co was known 
analytically, up to now (7|. Our result for C agrees with the one in Ref. f7j. The rest of the co- 
efficients have been calculated here explicitly for the first time. Our calculations put restrictions 
on the universal functions (see ||-||) entering in higher heat kernel coefficients. Finally, we 
must point out that, even if we only list a few number of them in Table 1, in fact our program 
allows for the calculation of any coefficient to any degree of approximation in the a dependence 
(notice that the coefficients are given by rational numbers, as it should be). However, to obtain 
each one of them as a closed function of a is outside the scope of our purposes here. 

Several of the first heat kernel coefficients are singular at a = 1. This is a manifiestation of 
the lack of strong ellipticity in the corresponding boundary value problem |], [<J . 
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A very interesting feature of our results is the vanishing of the a-dependent part of C3/2 (or 
a 2 ), while the heat kernel itself becomes singular at a certain (critical) value of the coupling. This 
means that the conformal anomaly and the one-loop counterterms are not modified by the pres- 
ence of the Chern-Simons boundary action, what could be a reflection of a general property of 
Chern-Simons theories. Since the fundamental explanation of very important physical phenom- 
ena of lower-dimensional QED, such as the fractional quantum Hall effect and high-temperature 
superconductivity, seem to rely very heavily on the fundamental structures provided by Chern- 
Simons theories, we conclude that the property we have here found (by looking to a particular 
model) could have important physical consequences for those subjects. This issue is presently 
under investigation. 
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Table 1: The first ten heat kernel coefficients. 



27rV 



-3/2 



8vr 3 / 2 (a 2 - l)- 1 

1 l-a 2 )- 1/2 (2 + a 2 -8a 4 )+ 117 



7r 2 r 



7T 



3/2 



2a 4 
59 



7T 

r 



45 
2 '3 



x tt 4096 
7t 3 / 2 / 117919 



32 



1339 11a 2 115 a 4 435 a 6 




1024 4096 16384 

2048 a 2 10112 a 4 816896 a 6 
+ + + 



+ 



15015 109395 14549535 

43053 a 4 14431a 6 



57455 6787 a 2 
+ 



393216 



371148101 



65536 
571904 a 2 



1048576 1048576 
1259008 a 4 310784 a 6 



+ 



7T 



116396280 2909907 37182145 42902475 
633 16417555 468537 a 2 64091a 4 515873 a 6 



160tt 201326592 4194304 16777216 33554432 



7T 3 / 2 7399265868279 997830656 a 2 386748416 a 4 



V 60235074900 5019589575 
68230270976a 6 



8562829275 



1933976154825 



+ 
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